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Abstract. 



The PZT material near the morphotropic phase boundary where a monoclinic 
phase was observed is studied in this paper. A theory of the tetragonal to mon- 
oclinic phase transition was developed within the Landau free energy approach. 
The order parameter is the electric polarization vector. The phase transition 
from the tetragonal to the monoclinic; phase is of the first order near the second 
order. The tetragonal and monoclinic phases are described as equilibrium states 
within this theory. Dielectric response near the phase transition from the tetrag- 
onal to the monoclinic phase was studied for the crystal and for the ceramic and 
polycrystalline material taking into account the diagonal 33 and the off diagonal 
31 components of susceptibility. The later diverges at temperature of stability 
boundary for the monoclinic phase, and it is zero in the tetragonal phase. Such 
a behavior is not observed in the effective dielectric static constant temperature 
dependence for ceramics and polycrystalline materials. We explain absence of 
a large contribution of the off diagonal susceptibility by a grain structure of 
the material and by percolation properties of this microcomposite. Above the 
percolation transition the off diagonal contribution to the effective susceptibil- 
ity is not present. Observed temperature behavior of the effective dielectric 
susceptibility is thus qualitatively explained. Effective susceptibility which is 
decreasing in tetragonal phase is in the monoclinic phase also decreasing, how- 
ever below phase transition temperature the decrease is slower. A quantitative 
theoretical behavior of dielectric susceptibility for both phases is found. Me- 
chanical response of PZT material to a uniaxial stress is studied theoreticaly for 
both phases. Observed temperature behavior of it, a well around the transition 
temperature from the tetragonal to monoclinic phase transition, is qualitatively 
explained. 
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1 Introduction. 



PZT materials are interesting ferroelectric materials studied for many years, [Tj 
to Near x ~ 0.48 a morphotropic phase boundary (MPB) exists It is 
a region between the tetragonal phase, which is on the Ti rich side, and the 
rhombohedral phase, which is on the Zr rich side. This region is associated 
with many interesting dielectric, piezoelectric and other properties of PZT j^j. 
The width of MPB varies with the homogeneity of the powders [71 to US] and 
the size of grains jJUl- Recently a new phase was found around MPB in the 
phase diagram. This phase is a monoclinic phase 5,_. Measurements of the 
elastic response ^T] on PZT have shown the presence of the transition from 
the cubic to tetragonal phase and a well at lower temperatures. The transitions 
from the cubic to the tetragonal and from the cubic to the rhombohedral phases 
are understood quite well, see in [Q. A transition from the tetragonal to the 
monoclinic phase which is not well understood will be studied in this paper. To 
describe the tetragonal to the monoclinic phase transition we need to identify 
the order parameter and then to use the Landau theory of phase transitions, 
this transition is of the first order near the second order. Dielectric response and 
mechanical response of PZT will be studied for the monoclinic phase and for 
this phase transition in this paper. We compare our results with experimental 
results, 11 and J^l- These materials are ceramics and polycrystals, they are 
composites and this fact does lead to a modification of the dielectric response as 
it is in the case of some ferroelasticdielectric microcomposites ^^l- Crystallites 
with different phases, tetragonal and monoclinic or rhombohedral, will not be 
considered here. The effect of different phases coexisting in PZT was considered 
in [H]. 

The paper has the following sections. In the following section we describe 
briefly PZT and the Morphotropic Phase Boundary problem. Then we study a 
phase transition from the tetragonal phase to the monoclinic phase. Transition 
temperature region between the monoclinic and the tetragonal phase is studied 
in more details. This enables us to study dielectric response near the phase 
transition. PZT material is usually ceramics or polycrystalline. Then its micro- 
composite properties influence the dielectric response which we study in the next 
section. Mechanical response of PZT near the phase transition from the tetrag- 
onal phase to the monoclinic phase is studied theoreticaly then. Dielectric and 
mechanical responses theoreticaly described in this paper for this temperature 
range qualitatively describe observed dependencies of them on temperature. In 
the last section there is a summary of obtained results in this paper. 

2 PZT and the Morphotropic Phase Boundary. 

The morphotropic phase boundary (MPB) in PZT is not well understood 0. 
There is evidence that this phase boundary is a two phase region, see in |14| and 
|15j . There are also measurements which show that this is a single phase region, 
or part of it is a two phase region (tetragonal and monoclinic) and part single 
phase region, see in ^1] where the author addresses the question of coexistence of 
different crystallographic phases within the same crystallite. Materials similar 
to PZT which have MPB are PMN PT and PZN PT [H]. It was found in 
[TH] that when the diffusion of Zr and Ti atoms in PZT is fast enough a two 
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phase region exists at this boundary. When the diffusion is not fast enough a 
nonequihbrium region forms and the classical phase diagram 12 appears in this 
material. 

To understand the nature of the MPB we will use the Landau theory for 
phase transitions. As it is noted in |14| the monoclinic phase as a intermedi- 
ate phase between the tetragonal and the rhombohedral phase may exists and 
that the tetragonal space group P4mm and the rhombohedral space group R3m 
both have as a subgroup the monoclinic space group Cm . Thus there exists |14) 
the possibility that a tetragonal to rhombohedral phase transition could take 
place via a monoclinic intermediate state. It is then very usefuU to understand 
properties of the monoclinic phase and of the phase transition from tetragonal 
to monoclinic phase. In ^Hl the authors have concluded based on local elec- 
tron diffractions, that the local structure of PZT is always monoclinic, and that 
the average tetragonal and rhombohedral structures are a result of shortrange 
to longrange ordered states. According to these authors these viewpoints may 
stimulate new constitutive models of ferroelectrics with compositions close to a 
MPB, and that so far, a constitutional model which deals with changes of the 
structural order is not available jl6| . It is then very useful to understand prop- 
erties of the monoclinic phase and of the phase transition from the tetragonal 
phase to the monoclinic phase. First principles calculations IT have shown the 
stability of the monoclinic phase in a narrow region of MPB when a random 
Zr/Ti cation distribution was taken into account. The stability of a monoclinic 
phase of space group Cm in between the rhombohedric and tetragonal phases 
was also explained using a LandauDevonshire approach |18| . The authors of 
[T^ predicted existence of two other monoclinic phases between the rhombohe- 
dral and orthorhombic phases, and the orthorhombic and tetragonal phase of 
perovskites, respectively. The later authors have found that the Landau theory 
of the sixth order in order parameter for the transition from the cubic to the 
tetragonal and to the rhobohedral phase cannot explain presence of a mono- 
clinic phase and thus they extend the theory to the eight order and then they 
are able to explain presence of a monoclinic phase. The monoclinic phase of 
PZT ceramics, Raman and phenomenological theory studies were done in 
The authors fl^) have found using Raman scattering the tetragonal to mon- 
oclinic phase transition in PZT ceramics near morphotropic phase boundary 
at low temperatures. They have found that the transition is characterized by 
changes in the frequency of lattice modes with the temperature. To discuss the 
stability of the monoclinic phase the authors jEl used the Landau Devonshire 
phenomenological theory. While these authors and to [23 were able to pre- 
dict the rhombohedric to tetragonal phase transition at approximately x=0.5, 
their calculations show that the monoclinic phase is not stable for any value of 
Ti concentration. Free energy as a function of Ti concentration for the rhom- 
bohedric, tetragonal and monoclinic phases was calculated by these authors at 
80 degrees of Celsia. According to ^Hl this result indicates that the free energy 
used in j^O] to is not suitable to describe the existence of monoclinic phase. 
To describe the monoclinic phase it is taken in (l9j into account an order pa- 
rameter, which is the monoclinic angle. It is noted by these authors that the 
monoclinic angle is not the true order parameter since it does not go to zero 
at the temperature of tetragonal to monoclinic phase transition, and they are 
using a modified order parameter in their Landau theory. They have observed 
that the value of the free energy of the monoclinic phase is now lower than that 
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of free energy of the rhombohcdric and tetragonal phases in a small region close 
to MPB and they noted that this result shows that the modified free energy 
correctly predicts the existence of a monoclinic phase in the MPB region, in 
accordance with the experimental observations. The vector of the electric po- 
larization is according to these authors not a correct order parameter to describe 
the phase transition from tetragonal to the monoclinic phase and they are using 
the third order parameter for the second order phase transition which enable 
them to describe this phase transition. In fact this third order parameter does 
not follow from the symmetry analysis on which the Landau theory of the sec- 
ond order phase transition and of the first order near the second order is based. 
A priori order parameter thus does not enable to understand the monoclinic 
phase and the transition from tetragonal to monoclinic phase correctly. We will 
study the monoclinic phase and the transition from tetragonal to monoclinic 
phase within the Landau theory. We note that the free energy expansion found 
in [201 to |2y for the cubic to tetragonal phase transition should be modified 
by terms dependent on the electric polarization and which are invariant in the 
tetragonal phase. The tetragonal to monoclinic phase then will be described 
correctly. It is driven by temperature and the primary order parameter is the 
electric polarization vector. The origin of terms dependent on the electric polar- 
ization is in interaction of polarization vector with fluctuations of other degrees 
of freedom which when integrated out in the tetragonal phase lead to invariants 
for this phase. This is consistent with the symmetry group analysis and group 
subgroup relations for possible second order and first order near to second or- 
der phase transitions, see in JJ, and In the author has studied the 
equations of compatibility, which have to be satisfied. He analyzed the three 
possible phase transitions: tetragonal rhombohedral, tetragonalmonoclinic, and 
monoclinicrhombohedral. He used the micromechanics approach. The author 
found, that the tetragonal rhombohedral transition may be compatible at a 
composition at the very Zrrich side of the MPB x k, 0.54. He notes that in this 
compositional region there are no experimental data for the tetragonal phase, 
so that it may not exist at this Zr concentrations. For the tetragonal to mono- 
clinic phase transition he has found that it is never compatible according to the 
accessible data and that the monoclinic to rhombohedral phase transition may 
be compatible at even higher Zr contents as tetragonal rhombohedral. He then 
discusses the incompatibility between the tetragonal and the monoclinic phase, 
and he has found that it is almost constant across the MPB, but relatively low 
compared to the other transitions. This may explain |14j . why tetragonal and 
monoclinic phases seem to coexist within the MPB region. 

We will study in our paper the tetragonal and monoclinic phases, the tran- 
sition between them we have found is of the first order near the second order. 
We will not study coexistence of these phases in our paper. 

The Landau expansion of the Gibbs free energy for the cubic phase was 
described in [211 to [21] and has the form: 

(1) 

+ani(Pf + P| + Pi) + aMPiiPt + ^2 ) + A'(^2 + ^3 ) + A'(^2 + ^3 ))+ 

+ai2^PlP2Pi] 
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here Pi (i=l,2,3) are the polarization components along the cubic axes, ai is 
the usual parameter from the Landau theory, (i,j — 1,2) and aijk (i,j,k=l,2,3) 
are the dielectric stiffness and highorder dielectric stiffness at constant stress, 
respectively. The coeffcient ai — aifl{T — Tq) is a temperature dependent 
coefficient. Here Tq is the transition temperature. The values of all coeffecients 
in 1^ and the solutions for the tetragonal and rhombohedric phases can be 
found in [201 to [23 ■ In the cited papers there is a Landau expansion Q and an 
expansion in other secondary order parameters. Contribution from the elastic 
tensor and its coupling to polarization is neglected in the expansion While 
the free energy expansion is for the cubic phase, in the case of the tetragonal 
phase ^3 the cations displacements lie close to the monoclinic [201] direction, 
which corresponds to the rhombohedral [111] axis. Thus the transition from the 
tetragonal to the monoclinic phase will be described by the free energy from 
(^. This phase transition is of the first order near the second order. However 
note that the free energy will contain contributions which correspond to the 
tetragonal phase symmetry now. One of such fourth order terms has the form: 



The term ^ breaks the cubic symmetry, thus it is not present in the corre- 
sponding free energy expansion for the cubic phase. The constant 6 is not 
known as concerning its value, we assume that it is not dependent on tempera- 
ture and pressure. However due to the fact that a transition from the tetragonal 
to the monoclinic phase is present in experiments we may assume that it has 
a very large value such that we obtain from the contribution (2) the equality 
Pi = P| at the ground state. This condition is used a priori for description of 
the monoclinic phase also in [2Z]. However its symmetry origin in [2Z| is not 
explained. The authors of [23 are taking into account a third order parameter 
in order to describe the monoclinic phase, which is the monoclinic angle. This 
approach is artificial. The case of 6 general in |(2l will be studied elsewhere as 
well as mechanical and dielectric response properties of the system for 5 general. 
Then we may find equations for the polarization components Pi and P3 only. 
Thus the order parameter space is two dimensional. 

Another of such terms has the form: 



The term breaks again the cubic symmetry, thus it is not present in 
the corresponding free energy expansion Q for the cubic phase. The constant 
T is not known as concerning its value, we assume that it is not dependent on 
temperature and pressure. Thus this term prefers the polarization to be oriented 
in the tetragonal axis. However we will not consider this term further in this 
paper. One may perform a systematic analysis of the free energy expansion 
corresponding to the symmetry of the tetragonal phase. We decided to use the 
terms and 121l in this paper because the coefficients of the terms in Q are 
well known and the term ^ gives zero contribution for the tetragonal phase 
equilibrium state P| > 0, = P2 = and for the monoclinic phase equilibrium 
state P| > 0, Pi — P2 > 0, and it has tetragonal symmetry. The case of 6 
and r general will be studied elsewhere as well as corresponding mechanical and 
dielectric response properties. 




(2) 




(3) 
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Decreasing temperature we may expect that the Landau theory based on l|T)l 
and (01 describes the transition from the tetragonal to monochnic phase. In the 
later phase we will have = P| | P3 . 

3 Phase Transition: Tetragonal Monoclinic Phase. 

The equations from which we will find the values of the order parameter in the 
equilibrium phase for a given temperature and pressure, i.e. of the polarization 
components = P| and P3 are: 

4aiPi + 4(2aii + ai2)Pi + (4) 

+4ai2PiP3 + 12( am + q;ii2)Pi + 
+4aii2PiP3* + 4(2aii2 + ai23)Pi'P| = 0. 

and: 

2aiP3+4aiiP|+ (5) 
+4ai2P3Pi^ + 6aiiiP3'^+ 
+8aii2Pf P| + 2(2an2 + ai23)Pi^P3 = 0. 

The solution of the equations Q and ^ for the cubic phase Pi = P2 = 
P3 = 0, for the tetragonal phase Pi = P2 = and P| | 0, for the rhombohedric 
phase Pi = P2 = P3 1 0, and orthorhombic phase Pi = P2 J and P3 = can be 
found. However these phases were already studied in [201 to 

In our model due to large constant S the polarization components fulfill 
Pi ~ P2 to for the equilibrium state of the monoclinic phase. Let us describe 
this monoclinic phase. To do this we introduce polar coordinates in the order 
parameter space (Pi,P3): 

Pi^Asm{(j)) (6) 
P3 = Acos{(j)) 

where A is the amplitude of the order parameter in this parameter space, 
and (j) is the angle of the order parameter in this parameter space. Note that 
for the order parameter angle equal to zero, (/> = 0, the polarization is oriented 
in the tetragonal direction. For the order parameter angle nonequal to zero, 
however smaller than ^, i.e. for ^ > > 0, the polarization is oriented in 
the monoclinic direction. For the angle 4> equal to j the polarization vector is 
oriented in the rhombohedric direction. For the order parameter angle equal to 
^ , the polarization vector is oriented in the orthorhombic direction. 

We have found from the equations Q to Q that the amplitude A depends 
on the phase angle cj) in the following way: 

= ^TT- (7) 

r + q sin [ip) 



Here r and q constants are defined as: 



and: 

2aii2 - ai23 

1 " ■ 

iaii - ai2 

We may introduce the new phase angle (3 by: 

TT 



(9) 



2 /? (10) 



and then the amplitude {T)) takes the form: 

= ^ (11) 

r + gcos2(/3) 

The angle /3 describes the tetragonal phase if it has value ^ . Near the phase 
transition from the tetragonal to the monoclinic phase the angle (3 < ^. From 
the equation lll|l we see that the amplitude A is almost independent on this 
angle. Near the phase transition from the tetragonal to the monoclinic phase 
we expect that the amplitude of the polarization vector does not change too 
much with temperature. This assumption corresponds with the experimental 
observations, to [SJ. Near the phase transition tetragonal to monoclinic 
phase the angle of the electric polarization changes mainly. We have found that 
at the transition temperature region the amplitude A is given by: 

^2 _ 3aii ~ ai2 ^^2) 

2aiii - q;ii2 + 0{(t))^ ' 

As we can see from H12|l the amplitude A of the order parameter is weakly 
dependent on the angle /3, or 0. However it is nonzero above and below the 
transition temperature i.e. we may assume that it is almost independent on 
temperature near the phase transition from the tetragonal to the monoclinic 
phase. Then we may use the constant amplitude approximation A = const, in 
our calculations near the phase transition temperature from tetragonal to the 
monoclinic phase. In this approximation we will assume that: 

^2 ^ _ 3an-«i2 (^3^ 
2q!iii - aii2 

is constant. This approximation is corresponding to breaking the cubic sym- 
metry for nonzero values of the angle 4>. Note that the free energy expansion 
l|T)l for the cubic phase is a part of the free energy expansion describing the 
tetragonal to the monoclinic phase, however the free energy will contain contri- 
butions which correspond to the tetragonal phase symmetry now. The lowest 
order term (O has now the form: 

SF = J dVi^A\cos{x)^ - sin(x)')' sui{^)\ (14) 

Here the angle x is the angle between the 1 and 2 components of the po- 
larization. For very large v constant this angle is 45 degrees in the monoclinic 
phase. We will assume that this angle has the value 45 degrees in the following. 
The discussion of a finite smaller value of the constant in (PJ will done else- 
where. From the equation and (jsj we obtain the following equation for the 
phase variable (ji: 

t\2ai -f Aai2A^ + 2(2aii2 -I- am)^^] + (15) 
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+i2[4ai + 4(aii + ai2)A^ + 8aii2A^] + 

+ [2ai + 4aiiA^ + 6aiiiA^] = 0. 

where t — tan(/3). The monoclinic phase is reahzed when the free energy of 
it wiU be lower than that of the tetragonal phase. Let us discuss now solutions 
of the equation H15|l . These solutions will correspond to extremes or inflection 
points of the free energy. The equilibrium state will be that solution which 
gives the lowest free energy. The equation is of the second order in the variable 
= tan^(/3). The discriminant D of this equation is: 

D = [4ai + 4(aii + aia)^^ + SaiiaA*]^ - (16) 

-4[2ai + 4ai2A^ + 2(2aii2 + ai23)A^][2ai + 4:aiiA'^ + 
+6aiiiA'^] 

and should be real. Let us consider the real case. In the case of imaginary 
discriminant D the monoclinic phase does not exists. Let us again consider 
the small amplitude A approximation. While this approximation is certainly 
correct near the cubic tetragonal phase transition, it is not clear whether it is 
correct near the tetragonal to monoclinic phase transition. However we will 
assume that the difference (2aii — 0:12) is small. The small constant amplitude 
approximation will be discussed numerically at the end of this section. In this 
approximation we may neglect the fourth order terms in A in the equation H16|l . 
In this case the discriminant D has the form of the square of difference of the 
coefficients of the zeroth order (let us denote it as c) and of the fourth order (let 
us denote it as a) in the equation H15|l . Then we have found that there exists a 
positive solution of the equation l|15|l for a > > c and for a < < c for . In 
other cases of a and c there is a tetragonal phase realized. The solution of 115|l 
in the cases a > > c and a < < c has the form: 



a 

or: 

n2 . C 



taniPf =1 - I (17) 
cot((/.)2 =1 ^ I (18) 



a 



which gives: 



or: 



tan /3f = \ 19 

^ot^"^) =1 ' o ' 42 (20) 

Here the amplitude A is constant in temperature in this approximation, it 
is given by H13(l . The equation l|19|l may be rewritten for small angles (3 in the 
form: 

=1 I (21) 

near the temperature Tm- The equation H20|l may be rewritten in the form: 
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near the temperature . Here the temperatures Tm and are defined 

as: 

Tm = To- (23) 

"1,0 



and: 



T^ =To — (24) 

"1,0 



The difference of both temperatures is: 

T^^_TN^_2i^i2^^n)A^ (25) 

"1,0 

which is a positive or a negative quantity. For the phase transition from the 
tetragonal to monoclinic phase this difference is expected to be negative. 

From the equation H22|l we see that the angle (f> is zero at temperature 
T^ . This temperature is identified with the highest possible temperature (the 
metastable) monoclinic phase may exist. Above this temperature there is the 
only solution of the LagrangeEuler equations, the solution in which the only 
nonzero component of the polarization, P3 , exists. This solution corresponds 
to the tetragonal phase. Below this temperature there exist nonzero components 
of the polarization, — P2 t P| . 

From the equation H22|l we see that the angle (p = ^ at temperature T — Tm 
. This temperature is identified with the lowest possible temperature of the 
metastable monoclinic phase to exist. Above this temperature there may exist 
nonzero components of the polarization, = P2 t P| . 

From the equation H22|l we see that at the angle </> = ^ at temperature Tm = 
Tm+tn ^YiQ rhombohedral phase angle appears. Components P^ — Pi = Pi 
are equal as concerning their amplitude at this temperature. 

The temperature Tn at which the rhombohedral phase angle value becomes 
present is given by: 

^^_^^^Jau+an)^^ (26) 
"1,0 

The calculations above were done for the approximation in which we ne- 
glected the terms of the order 0(4) in the amplitude A. The angle devel- 
opment was calculated for temperature interval T]\f < T < T^ , i.e. for the 
interval from the highest possible temperature T^ for which the monoclinic 
phase may exist, to the lowest temperature for which the monoclinic phase may 
exists. The temperature at which the rhombohedral angle value of realizes is 
at the middle of this interval. 

We have calculated the numerical values of the amplitude for the con- 
centrations 50/50, 40/60, 30/70 and 20/70 using the numerical values of the 
constants in the Landau expansion for the free energy from [20] to [53] . We 
have found that it has a negative value 3.690 in the corresponding units for the 
50/50 ratio, and the positive value -1-0.834 for the 40/60 ratio increasing further 
its value to -f-6.039 and +11.880 for 30/70 and 20/80 ratios respectively. For 
the 10/90 ratio it has a negative value 4.507. We have found that the constant 
small amplitude approximation near the boundary 50/50 is a good approxima- 
tion for PZT and that is a small positive parameter in this region where it is 
changing its sign. The amplitude A has nonzero value at tetragonal phase. The 
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phase transition from tetragonal to monoclinic phase is of the first order near 
to the second order. Thus the negative numerical value of the amplitude for 
the monoclinic phase for the concentrations 50/50 means that the monoclinic 
phase does not realize for this concentration. Note that the amplitude is 
calculated for the monoclinic phase. Decreasing temperature the free energy of 
the monoclinic phase is expected to increase, and at some point it may become 
larger than that of the rhombohedral phase. 



4 Transition Temperature Region between the 
Monoclinic and the Tetragonal Phase. 

Let us consider the free energy of the tetragonal and the monoclinic phase. 
Comparing them for the equilibrium states we obtain the transition temperature 
between the monoclinic and the tetragonal phase. The later phase is that phase 
for which = P| = and P| | 0. The former phase is that phase for which 
= P| I and P| = 0, the components 1 and 2 of the polarization vector 
are different from the component 3. The transition temperature between the 
tetragonal phase and the cubic phase, in which the polarization is zero, is: 

2 

Tt-c = To+ . (27) 
4Q;i,oaiii 

As we see this temperature is different from temperature Tq which contains 
the coefficient ai . The phase transition is of the second order near the composi- 
tion 50/50 and of the first order near the second order for higher concentrations, 
see 120] to 123. One can use the concentration dependence of the coefficients a 
in the free energy expansion and find how this transition temperature, and also 
other quantities, is changing with the concentration, [201 to [21] ■ This is not 
our aim here. The free energy of the monoclinic phase will be compared with 
that of the tetragonal phase. The free energy of the monoclinic phase is given 
within our approximation of v large by: 

'-0 



I dV[a^A'j^+ (28) 



Here io = tan((/)). The free energy of the tetragonal phase has within our 
approximation the form: 

F ^ j dV[aiA^ + aiiA^ + aixiA% (29) 

In our approximation we have in the tetragonal phase P3 = A^ . The 
amplitude A^ in the tetragonal phase is given for PZT in the range near MPB 
by: 

^T2 ^ -2q:ii + ^{2aiiY - 12o;iiiai ^^^^ 
6aiii 
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The transition temperature from the tetragonal to the monochnic phase is 
found from free energies of the tetragonal phase and of the monoclinic phase 
comparing them. At this temperature they are equal. Above this temperature 
the free energy of the tetragonal phase is lower than that of the monoclinic 
phase, below this temperature the free energy of the tetragonal phase is higher 
than that of the monoclinic phase. For the monoclinic phase we will use the 
small constant amplitude approximation as described above at this region of 
temperatures. Then we obtain for the transition temperature the equation: 

SD-^=4al,-2ai.A-^'-^ (31) 
ail ail 

where the discriminant D is given by: 

D = (2aii)2 - Uaiiiai (32) 

Note that in between 40/60 and 50/50 concentrations there is present a 
large change of the constant y = ^ 4^"^^)"^ where in ai we take the temperature 
for the transition from the tetragonal to monoclinic phase as 358K. For the 
concentration 40/60 this constant has its value 1.179 and for the concentration 
50/50 the value 60.387 . The discriminant D above can be written as (2q;ii)^(1 + 
2y). Nearby the MPB the constant y is expected to be large. Then we obtain 
from (32) that the transition temperature Tt-m from the tetragonal to the 
monoclinic phase is given by: 



J^T-M = -Lq - 7, (l + Q!i2A . ^ — ) (33) 

3ai,oaiii 4aJi 



To estimate the numerical value for this transition temperature we need the 
values of the constants in (|33|l for the concentrations near MPB for concentra- 
tions near 48/52. These are not available. For 50/50 concentrations we obtained 
from that temperature Tt-m in is negative. It means that at this con- 
centration there is no monoclinic phase. This corresponds with experimental 
findings up to now. For concentrations 40/60 our approximation of a small con- 
stant amplitude is not valid. Moreover the rhombohedric phase is preferred at 
this concentration as that phase which evolves from the cubic phase decreasing 
temperature. 



5 Dielectric Response Near the Phase Transi- 
tion. 

The dielectric response of PZT in the tetragonal and monoclinic phases will be 
calculated from the free energy given by and Q in the approximation of large 
ly applying an external electric field E in the direction 3. The external field will 
interact with the polarization, and this interaction gives a contribution — E.P 
to the free energy density. One can ask how PZT responses to the field which 
is parallel and perpendicular to the tetragonal phase polarization direction in 
the monoclinic phase. One can also ask how PZT responses to the electric field 
in the monoclinic phase where there exist polarization components (1,2) which 
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are absent in the tetragonal phase. Dielectric response will be calculated from 
the total free energy for the tetragonal phase. 

Let us first consider how PZT responses to the fields which are parallel to the 
tetragonal phase polarization axis in the tetragonal phase and in the monoclinic 
phase. Let us assume that there is a nonzero electric field in the direction 
(001) in pseudo cubic notation. 

In the tetragonal phase the field parallel to the polarization induces changes 
in the amplitude of the polarization vector, the monoclinic angle (j) is zero and 
remains zero for this field. The change SA in the amplitude A induced by this 
field is: 

From the equation (ref34) we see that there is present a diverging behaviour 
from the cubic to tetragonal phase in the susceptibility 633. Note that D is 
defined in the equation (|32|l and the amplitude A in the tetragonal phase is 
given by (|30|) . This field induces a change of the polarization in the direction 
(001) and does not induce a change of the polarization in the direction (100) or 
(010). 

Near the transition from tetragonal to monoclinic phase the amplitude A in 
the monoclinic phase is moreless constant, we will use this approximation. 

The angle (j) is nonzero in the monoclinic phase, thus we may expect that 
nonzero electric field E3 induces a change not only in the direction (001) but 
also in the directions (100) and (010). From the equations l@J and © we find 
that the polarization components are related by the equation: 



1 + r{P,' + P|) + qPi = ^p2'^_%^p^ - (35) 



In the polar order parameter space A, (j) we obtain from H35|l: 

l + rA^ + qA^ sin^ (0) = — ^^fi . (36) 

^3COS(0)COS2(20) 

As we can see from the equation (|36ll the amplitude change induced by the 
external field E3 is coupled to the angle (j) change induced by this field. 

We have found that the dielectric susceptibility £33 in the monoclinic phase 
is given in the monoclinic phase within our approximations by a contribution 
due to the induced change of the angle: 



Ph ni+tl?) 
^33 



q,/l]Tr{2-7tl + 4) 



(37) 



and by a contribution due to the induced change of the amplitude: 

-2r 



■^33 



3^/1 + 4 



(38) 



The total dielectric susceptibility e^l for the monoclinic phase is given within 
the constant amplitude approximation by the sum of (|37|l and 138(1 . 
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We have found that the dielectric susceptibiUty 613 in the nionoclinic phase 
is given within the constant ampUtudc approximation in the monochnic phase 
by: 

£13 = £33 + (39) 



Here we denoted to = tan((/)) and e|" 



£^3 = (40) 



■^33 ■ 

2r 

y 

It is easy to find the susceptibiUty of the tetragonal phase as we described 
it above, see also El to 123 It is temperature dependent and diverging at the 
trantision temperature from the cubic to the tetragonal phase for the second 
order phase transition, and almost diverging for the first order phase transition 
near the second order. 

Near the tetragonal to monoclinic phase transition temperature we have 
found that dielectric susceptibility ei3 is zero in the tetragonal phase: 

£13 = 0. (41) 

Near tetragonal to monoclinic phase transition temperature we have found 
in the monoclinic phase the dielectric susceptibility £13 to be given by: 

ei3 = £^3[^^^ ] (42) 




The dielectric susceptibility £13 diverges at the temperature . 

As we can see the susceptibility £33 in the direction of the applied electric 
field is nonzero in the tetragonal phase and in the monoclinic phase. Decreasing 
temperature it is decreasing in the monoclinic phase in a slower way than in 
the tetragonal phase. The slope of decreasing susceptibility in tetragonal phase 
changes at this transition temperature to smaller one. This behavior is seen in 
the experiment, see The susceptibility £33 has its value at the temperature 
Tn ■■ 

£33 = £i3[^ + ^=] (43) 



which leads to the value ^ — 0.906 at this temperature. As we can also see 

^33 

the susceptibility £13 is zero in the tetragonal phase and in the monoclinic phase 
it is either increasing either decreasing function in temperature depending on 
the sign in the equation H39|l . For the 40/60 concentration ratio we have found 
that it is a negative quantity. Applying the field in the monoclinic phase 
there is a contribution to the susceptibility due to the change of the angle 
<f> in such a way that the polarization vector is tilted to the tetragonal axis 
(001), thus the £13 susceptibility is expected to be negative. Note that the 
constant r is independent on temperature in our approximation, however it is 
dependent on the PZT solid solution concentrations ratio. So there is a jump in 
this component of the susceptibility. There is also a jump in the susceptibility 
component £13 present for the crystal case studied. 
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6 PZT Ceramics Their Micro composite Prop- 
erties and Dielectric Response. 

The ceramic PZT is composited from crystallites which have polar axis oriented 
in different directions and which have different shape, and from a free volume. 
There exists a possibility of coexistence of crystallites with different phases 
tetragonal and rhombohedral at and nearby MPB which leads us to the question 
whether such coexistence does not lead to a possible enhancement or decrease 
of the dielectric response "SF in this material. The effect of different phases in 
different crystallites in PZT was considered in [ij. According to this work we 
can consider PZT as a microcomposite. Let us calculate the effective response 
of the PZT microcomposite. It consists of grains which have the tetragonal axis 
(001) oriented in a different ways. The effective dielectric susceptibility of the 
microcomposite eg// is found within the Effective Medium Approximation j28| 
from: 

.^^//-^i , (l_^)i£Z^^^O. (44) 



Here the susceptibilities ei and e2 are susceptibilities for two types of ma- 
terials (grains) in the microcomposite. The concentration x is a concentration 
of the first type of the material, the concentration Ix is a concentration of the 
second type of the material. Let us assume that the first type of material is that 
part of grains which has its tetragonal axis oriented in the direction of the field. 
Let us further assume that the second type of material is that part of grains 
which has its tetragonal axis oriented perpendicular to the direction of the field. 
This is a simplified model of the PZT ceramics, we expect that it gives a qual- 
itative picture of the microcomposite response. The effective eg// constant is 
given in the tetragonal phase, where the off diagonal dielectric susceptibility is 
zero, given by the £33 constant. In the monoclinic phase the offdiagonal dielec- 
tric constant is large, diverging at . It can be easily found from (|44|l that 
the effective dielectric constant for concentrations x of the first type of material 
above the critical concentration for percolation of spheroidal grains Xc = | is 
given by: 

= ra- («' 

Thus it is a finite response, determined by the first type of grains. Note that 
the closer we are to the critical concentration the larger is this effective dielectric 
constant. It is this behavior which corresponds to the observed behavior, see in 
|12j . Below the critical concentration the effective: 

'^efS = ^{X-Xc). (46) 

This effective dielectric constant is large, almost diverging depending on 
divergence of the dielectric constant £13. 

In the tetragonal phase we have found that the effective susceptibility is 
given by: 

£e// = £33- (47) 
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7 Mechanical Response Near the Phase Transi- 
tion. 



Let us consider a crystal of PZT which is under the uniaxial stress a. We will 
use symmetry considerations of Haun [21] to [21] to find a Gibbs free energy 
expansion. However in this case we need to consider the coupling between the 
elastic tensor components (strain) and the polarization vector components. The 
nonzero uniaxial stress interacts with the M3 component of the elastic tensor. 
We are using notation 1 to 6 for tensor components of the strain tensor u. This 
component interacts with other components of the elastic tensor (strain) and 
with the polarization components in a way described by the part of the Gibbs 
free energy which has the same form as the part used by Haun et al. (20j to 
for the Gibbs free energy with stress aij, i,j=l, 2, 3. The thermodynamic 
potential 5F (the Gibbs free energy) expansion in strain is used in the same 
way as in Landau and Lifshitz book |25| : 

^ j dV[]^Cix{ul + u\ + u\) + Ci2(wiU2 + + U3"i) + (48) 

+ ic44(w4 + "5 + ^l) + + + -"3^3 ) + 

+gi2(Mi(F| + Pi) + ^2(^3' + Pi) + "3(A' + Pl))+ 

+qM{u4,P2P3 + U5P3P1 + UePlP2) - -"30-3]. 

The uniaxial stress in the direction 3, (T3 = tr in the following. The elastic 
tensor components (strain) have their values from the following equations: 

C11U3 + Cuiui + U2) + quPi + qi2{P? + Pi) - - 0, (49) 

and: 

CiiUi + Ci2(m2 + us) + qiiPl + qi2{P2 + Pi) = 0, (50) 

and: 

C11U2 + Ci2(u3 + "1) + 911^2 + 912(^^3 + ^i") = 0. (51) 
When the coupling between the polarization and the strain is absent, i.e. 
when qi2 = qn = 0, then the deformations in the 1 and 2 directions have 
the value, ui = 7/2 = — M3 q'^^q^^ , and in the direction 3 it has the value 
U3 ~ — , When the coupling is present, qi2 = ^0 and qn | 0, then we 

obtain for the component u^: 

U3 = - quPi - qi2 {P'l + Pi)+ (52) 

+ r ^I'r + ^2 ) + qi2iPi + Pi + 2Pi))]. 

Oil + 012 

As we see from ()52|l the strain tensor component U3 depends on the stress 
(T and on the components of the polarization vector. In the cubic phase Pi — 
P2 = P3 = and the tensor component U3 is given by: 

U3 = ^-^2 cr- (53) 
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In the tetragonal phase Pi = P2 ^ and P| > the elastic tensor compo- 
nent U3 is given by: 

1 9/^ 
"3 = — -W - (911 - 7^-^)Pi]- (54) 

In the monoclinic phase Pi = P| J P| = the elastic tensor component U3 
has the form in the A, <f> space: 

"-^in^ Cii(Cii+Ci2)'' " 

+ 7T^57^('ZiiA^t|t, + 'Zi2A^)]. 

'-'12 + L'll -L + ^0 

In the tetragonal and in the monoclinic phase we have the change of the 
strain tensor Su^ induced by the (small) uniaxial stress ti, for zero coupling 
between the strain and the polarization: 

<5w3 = ^~^cP 

*^ll(^ ~ Cii(Cii+Ci2)) 

In the tetragonal and in the monoclinic phase we have the change of the 
strain tensor components 5ui and 5u2 induced by the uniaxial stress a for zero 
coupling between the strain and the polarization: 

5ui=5u2 = ^^7^ — . (57) 

C,, J ii:i2_ Oil + C12 

In the previous calculations of the strain tensor we have neglected the influ- 
ence of the coupling between the strain and the polarization and we obtained 
the relations 1)56(1 and ((57(1 . We see from them that the strain tensor compo- 
nents are temperature independent. However this coupling has to be taken into 
account, in experiments 1%, we see that the strain is temperature dependent. 
From the equations 156(1 and 1(571) we define X33, 0:13, and 0:23 by U3 = a;33tT, 
ui — xi^a and U2 — X23cr. Here 0:13 — X23 ■ Taking into account the coupling 
between the strain and the polarization, i.e. the nonzero coefficients gn, 912 
and 544, the ai coefficient renormalizes to ai i and ai.s coefficients, where: 

Q!i,i = ai + (qiixis + gi2(a;i3 -|- 0:33)) (58) 

and where: 

ai,3 = ai + (9112:33 + '7i2(a;i3 + 2^23)) (59) 

In the monoclinic phase, taking into account the coupling between the elec- 
tric polarization vector and the elastic tensor (strain) components, we obtain 
the following Lagrange Euler equations for the components Pi and P3 : 

4ai,i -t- 4(2aii 4- ai2)Pf + 4ai2P3^ + 12(aiii + aii2)Pi^ + (60) 

-h4aii2P3^ + 4(aii2 + ai23)Pf P| = 
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and: 

2ai^3 + 4aiiP| + 4ai2Pf + GamPi + (61) 
8an2)P?Pi + 2(2aii2 + am)^* = 

From the equations H6()|l and (|61(l we obtain a relation between the stress 
and the polarization vector, neglecting terms A'^ in the small amplitude approx- 
imation, expressed in the polar coordinates A, cj): 

Q(T-4(2Q;ii+Q;i2)^^COs((?i) =0, (62) 

where the coefficient Q is defined by: 

Q = 4(giia;i3 + ^12(2^33 + 2:23)) - 4(1711X33 + Qi2(x23 + X13)). (63) 

As wee see from the equations (162(1 and ((6311 the stress is related to the am- 
plitude A and to the angle (f) when the coefficients qn and 512 are nonzero. The 
stress induces a change in the angle (f) which leads to a temperature independent 
renormalization of the component of the inverse compliance in our approx- 
imation. For the constant amplitude approximation in ((63(1 we obtain for the 

component of the inverse compliance in this case within our approximation: 
-f- - ^ — [1 - (2^12 - gii) ,, , + . (64) 

Cii(Cii+Ci2) 



2Ci 



2 



-qii) 



Cii + ^12 8(q;i2 — 2q!ii) 

From the equations ((62(1 and 1(63(1 it follows however that the stress induces 
a change in the amplitude A such that it leads to a temperature dependent 
renormalization of the component of the inverse compliance in our approx- 
imation. For the constant angle <j) in 1(63(1 we obtain for the component of 
the inverse compliance in our approximation: 

faf ^ J 1 . _ (gii + qi24 + c;^(gii)^§ + gi2) 

S<T Ciil^2^ / 4(2aii-ai2)(l-t2) 

(65) 

Thus a temperature dependent strain component in the monoclinic phase is 
found due to a change in the amplitude of the polarization vector and due to 
the coupling between the strain and the polarization. 

We can write the quantity in a general form appropriate for fitting 
experimental curves: 



-f— = 7=^ I + T0 + — -2-]- 66) 

OCT Gil 1 _ 2 1^ 1 ^ ^0 

C'ii{Cii+Ci2) 

We see from the equation (|()^ that the quantity is finite as the tem- 
perature T increases to the temperature T-^ , the temperature at which the 
monoclinic phase becomes unstable for higher temperatures. With decreasing 
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temperature the inverse of the function, which is measured in experiments 
and |12j . becomes larger in the monocHnic phase as can be seen from (Ifi5ll . 

In the tetragonal phase the inverse of the -j^, which is connected to the 
Young modulus, is calculated from its free energy expansion taking into account 
the strain contribution, the strain polarization component coupling and the 
stress strain contribution. We obtained that the strain component U3 has the 
form for nonzero uniaxial stress: 

U3 = ^-TTT^ [cr- (67) 

^11 Cii(Cii+Ci2) 



CJ 1 / 

-{qii - qi2-p, fr^a + J (4aii)^ " 12aiii(ai + crpi)))]. 

L-ii + <-'i2 oaiii V 

Here = an - gii(gii - gi2 c^^+k^ ) + 'JiaCcI^pfe hur—){qii 



and pi - i_(gn-9i2cnife)+9ii + 9i2(-c7^— 

With increasing temperature T the inverse of the see (I67|l . which is 

measured in experiments, increases in tetragonal phase. 

Temperature dependence of the inverse of the —g^ below the transition tem- 
perature from monoclinic to tetragonal phase is such that it increases as temper- 
ature T decreases. The inverse of the quantity is increasing with increasing 
temperature above the transition temperature from monoclinic to tetragonal 
phase. This behaviour is expected to occur only in the neighbourhood of the 
Morphotropic Phase Boundary, where we expect that (7ii-l-2gi2 is positive. Here 
qii w Qii and gi2 ~ Q12, the quantity Qn and the quantity Q12 are depen- 
dent on the composition of PZT, see in [221 • This behaviour of the mechanical 

response function corresponds qualitatively to the known experimentaly 



measured dependences, see in dl, and J2i- 



8 Summary. 

A phase transition between the tetragonal phase and the monoclinic phase was 
described within the Landau theory of phase transitions of the second order 
and of the first order near the second order. We have found that the transition 
from the tetragonal to monoclinic phase is driven by the electric polarization 
vector. So it is not necessary to introduce an artificial order parameter for this 
phase transition. The transition is of the first order near the second order. The 
tetragonal phase and the monoclinic phase equilibrium states were found from 
the free energy expansion for the tetragonal phase symmetry. We have used the 
constant amplitude approximation to describe the monoclinic phase and the 
phase transition from tetragonal to monoclinic phase. Such an approximation 
corresponds to the fact that the polarization vector does not change its ampli- 
tude too much during this phase transition, which was observed in [Tg. The 
transition is studied in the small constant amplitude approximation because a 
numerical estimate of the amplitude using the numerical values data for PZT 
free energy expansion coefficients in the cubic phase, as found by Haun et al. 
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|2U| to |24j . has shown that it is a small quantity in dimcnsionlcss units. To 
estimate the numerical value for the transition temperature we used the values 
of the constants in l|33() for the concentrations near MPB for concentrations near 
48/52. For 50/50 concentrations we obtained that temperature Tt~m in l|33() is 
negative, at this concentration there is no monoclinic phase. This corresponds 
with experimental findings up to now. 

Applying an external electric field in the tetragonal direction and an external 
uniaxial stress in this direction we studied in this paper dielectric and mechani- 
cal responses of the PZT crystal in the tetragonal and in the monoclinic phases. 
To explain observed dependencies in JI] and |E1 of the static effective dielec- 
tric susceptibility we had to take into account that samples on which dielectric 
response was measured were either polycrystals or microcomposites. A model 
for such a microcomposite was formulated in We have found that large 

off diagonal dielectric response present in the monoclinic phase near the transi- 
tion from the tetragonal phase due to its divergence at the temperature of the 
stability of the monoclinic phase is not seen in the experiments because of the 
polycrystalline and of microcomposite nature of samples in which the concentra- 
tion of the components with off diagonal dielectric susceptibility component are 
below the percolation threshold. Note that according to our theory one should 
observe at least a kink in this component in a pure crystal. We have succeeded 
to explain qualitatively the observed temperature behavior of the effective di- 
electric response, the formulas describing quantitative this behavior have been 
found in this paper, however a fit of them to experimentally observed curves 
will be done elsewhere We have studied the effective dielectric response of 
PZT polycrystals and ceramics using the Effective Medium Approximation. 

The mechanical response was also studied. The mechanical response ten- 
sor 33 and 13 (in the shortened notation 16) was found for a crystal for the 
tetragonal phase and the monoclinic phase. Its temperature behaviour corre- 
sponds to that observed in the experiment, there is a well around the transition 
temperature. 
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